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The problems of implementing the JAWS wind shear data are discussed in this pre- 
sentation. First of all, I will describe the data sets from the point of view of 
utilizing them in an aircraft performance computer program. Then I will describe 
some of the problems of non-standard procedures in terms of programming the equations 
of aircraft motion when the effects of temporal and spatially variable winds are 
included. Finally, I will show you some of the computed effects of the various wind 
shear terms. 

The specific tasks to be performed under NCAR contract are listed in Figure 1. 

The collection and processing of the dual and triple Doppler returns have been des- 
cribed in references 1 and 2. The processing of the Doppler radar signals resolves the 
data into a rectangular array of three Cartesian velocity components. The array is 
updated every two minutes, resulting in the fourth, dimension, time. Thus we have, 
essentially, blocks of data that represent every two minutes of the microburst 
phenomenon taking place. 

Figure 2 illustrates the origin of the coordinate system used on the NCAR data 
tapes. The origin of the x,y,z coordinate system is located at CP-2. The x coordi- 
nate is measured positive toward the east; and the y coordinate is measured positive 
toward the north. The z coordinate Is positive in the vertical direction upward. 

When implementing the wind data into the airplane equations of motion, one must be 
aware that the data are provided in an earth-fixed coordinate system. 

At present, we have looked at three data sets. The primary data set we have 
looked at is the August 5 data set. It was measured in the region indicated in 
Figure 2. The numbers in parentheses are the x, y coordinates of the region in 
kilometers relative to CP-2. The June 29 data set also shown in Figure 2 is another 
data set that came out earlier. These data are not processed in exactly the same 
manner as the August 5 data and they may show more severe wind shear effects than are 
real. We want to caution you to not use this data set because all the smoothing tech- 
niques and mass balance verification procedures currently perfected by NCAR have not 
been applied to the June 29 data set. The July 14 data set was measured at the loca- 
tion indicated on Figure 2. A great deal of work with this data set has not yet been 

carried out. Never before have we measured velocity to the resolution achieved by 
JAWS in a volume of space that is, essentially, 12 km square and 2 km high. We have 
wind velocity entirely throughout this volume element. The data is provided to us on 
data tapes with an established grid system. The grid spacing for the August 5 case is 
150 m in the horizontal direction (both lateral and longitudinal) and 250 m in the 
vertical direction (see Figure 3). 

In the 12 km square by 2 km high volume element, there are 81 grid points on each 

side, plus 9 grid points in the vertical direction. One of the problems that scares 

everybody immediately when you look at this data set is that you have 81 by 81 by 9 
grid points; and if you have three wind components at each grid point, you end up with 
approximately 177,000 data points. Storage capacity of a computer may thus begin to be 
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• Prepare from the JAWS data, four-dimensional comouter 
models on microburst velocities as input for flight 
simulator models. 

• Incorporate the new four-dimensional wind shear models 
into numerical computations to determine critical wind 
shear severity thresholds, access the scales of motion 
which lead to dangerous aircraft response, determine 
the relative importance of the horizontal versus the 
vertical wind speed component, define the test flight 
deterioration parameters, and evaluate operational pro- 
cedures for use in wind shear encounters. 

• Review flight simulator theories and aircraft equations 
of motion to assure compatibility of the wind shear 
models with existing simulator capabilities and compu- 
ter storage capacity. 

• Combine the instrumented aircraft high-frequency wind 
speed measurements with Doppler radar data to provide 
a meaningful turbulence model for addition to the wind 
shear simulation input models. 


Figure 1. Objectives of current research 
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Figure 2. Location of Doppler CP-2, CP-3, and CP-4. The origin 

of the XYZ coordinate is located at CP-2. The location 
of the August 5, 1982, and June 29, 1982, data sets are 
indicated by the shaded areas. 
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a problem. However, I think this is not true for the NASA simulators. In turn, at 
FWG Associates we have had problems with our very small super micro-computer (40 
kilobytes of storage) in storing that much data. Storage capacity on modern computers 
is not difficult to obtain. In the event that storage becomes a problem, we have 
selected volume elements of the data set, which I will describe to you in a moment. 

Recall that we are also talking about turbulence parameters for each grid station, 
such as turbulence intensities and length scales, if we can extract these from the 
radar signals. In turn, if you are interested in heavy rain, we can provide radar 
reflectivity. Thus, there may be as many as six parameters for each grid point. The 
storage capacity is thus doubled. The question of whether we can handle storage 
capacities of these magnitudes is one of the issues we want to discuss. Another 
problem which was identified was not so much storage as interpolation to find the wind 
speeds at the aircraft position at 64 times per second. This, I believe, is the typ- 
ical rate for most simulators. The procedure we currently use is a volume weighting 
procedure as illustrated in Figure 3. If the airplane is at point p, we simply weight 
velocity at a neighboring grid point with the diagonally opposite sub-volume element. 
Thus, we interpolate the winds at the point and when the airplane moves to the next 
point, we repeat the interpolation. This would seem to take considerable time on a 
computer, but actually this represents only about eight or nine lines of computer 
programming, and is not difficult to do (see ref. 3). Moreover, NASA is doing it in 
real time on their simulator. Thus, I believe, handling this large amount of data is 
not difficult, and the urgent need for a simple analytical model everyone is talking 
about is not that pressing. I firmly believe that simply storing the full volume data 
set is the simplest model you can get. By the time you come up with a model that is 
representative of true microburst wind shear, it is just as easy to use realistic data 
that have been measured in the field. 

Now let me address another issue which I hope will be discussed in detail. I feel 
it is very important to put the spatial derivatives or wind gradient into the aircraft 
equations of motion. There are nine spatial gradients of winds (Figure 4). These have 
never really been considered in previous analyses of aircraft motion, partly because we 
did not have the information available. I will show you later that the gradients do 
enter the analysis, and in some cases, we think can be very significant. Therefore, 
in addition to interpolating for the three wind speed components, you must also inter- 
polate for the wind gradients when using JAWS data sets. Again, this interpolation is 
not excessively time consuming (ref. 3). 

In order to assist you, the user of the JAWS data sets, we have identified paths 
through the full volume set along which wind shears occur. We have classified these 
as to severe, moderate, and light wind shear. The paths were initially selected by 
inspection. Figure 5 is a plot of the horizontal winds at ground level. The vectors 
indicate the direction of the wind and the size of the vector indicates the magnitude. 
The center of the microburst is clearly visible. We picked flight paths where a 
strong head wind changed to a strong tail wind along the path. If, for example, you 
fly through this data base, say along YZ, you have a strong head wind shifting to a 
tail wind. We picked a number of paths just by inspection. The path IJ was selected 
because it is interesting in that a relatively strong lateral wind shear occurs. 

Also, we selected some paths such as GH that we knew were relatively benign, but 
still challenging. We analyzed the particular wind field by conducting a computer 
simulation of aircraft performance along the respective paths. The aircraft was 
trimmed for a 3° approach along this flight path. Simple control law algorithms 
were used to maintain the approach path. The runway can be moved relative to the 
wind, or the wind relative to the runway, whichever way you like to look at it. A 
3° approach along flight path AB, for example, with the runway at any position 
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Figure 4. Wind vector gradients appear in the 
governing equations of motion for 
the aircraft (nine derivatives) 
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relative to where the center of the microburst is encountered can be simulated. We 
also simulated aircraft takeoff, with the position the aircraft flies through the 
microburst center again being adjusted. We ran several computer analyses and 
picked out what we classified as severe, moderate, and weak wind shear cases. The 
severe case on approach is one that our computer simulation, using relative simple 
control laws, was unable to fly; in other words, those cases in which the aircraft 
encountered the center of the microburst at a low altitude (roughly 100 m above the 
ground). In these cases, regardless of the fact that the control algorithm commands 
15° to 20° pitch and full thrust to go-around, the aircraft could not fly and 
crashed; we called these cases moderate. We called those paths weak where a 
go-around was possible without too much difficulty, but where there was enough wind 
shear to make the control laws move around significantly. 


A description of the paths and the nomenclature used to define them is shown in 
Figure 6. We trimmed the aircraft at roughly 2,000 ft. outside the data set. Then 
we picked some coordinates for each path which would represent the point at which 
the center of the microburst occurred. We measured those coordinates relative to 
the northwest corner of the data set. The coordinates of the NW corner shown on 
the figure are relative to CP-2. A new coordinate system for this specific data 
set is defined such that x 0 is the distance toward the east to the center of the 
microburst; y 0 is the distance toward the south to the center of the microburst, and 
Zq is the height at which the airplane would pass through that microburst if it were 
able to maintain the 3° glide slope along the designated path. For example, taking 
line AB shown in Figure 5, the runway is oriented along this path and the aircraft 
is trimmed to make a 3° glide slope; x 0 and y 0 then indicate the position of the 
microburst center relative to the end of the runway, and z 0 , which is a very 
critical parameter, is the height at which the airplane would pass through the 
microburst. The angle shown on Figure 6 is the angle of the path relative to 
the x n axis. 


To reduce the magnitude of the data sets for the convenience of the user, 
we picked particluar flight paths, for example, path AB in Figure 5, and 
prepared sub-volume data sets called 1-, 2-, and 3-plane models. The sub-volume 
data sets were constructed by passing a vertical plane and two parallel planes 
on each side separated by 500 ft. through AB. Data from the full volume data 
set were then transformed to these data planes. The longitudinal velocity, 

Wy , relative to the plane is the wind speed component along the direction of 
the particular flight path (i.e., AB in this case), is the lateral wind which 
is perpendicular to x, and w z is measured in the positive z direction upward. 

The center of the microburst (or roughly the center) is the origin of the 
coordinate system. Figure 7 schematically illustrates how the data are tabu- 
lated on a corridor data set. 


Figure 8 is a portion of a corridor data table for a 3-plane data set. A 
corridor along path AB is shown. The wind shear along this path is classified 
as severe, the origin is at the center of the microburst, the planes are separ- 
ated by 500 ft. Plane 1 is 500 ft. to the side of the center plane and 3 is 
500 ft. to the other side. Plane 2 runs directly down path AB (see Figure 5). 

The first column is the position along the path measured from the center of 
the microburst. The next three columns are the wind speed components in the 
x,y, and z direction, respectively. Notice that x is negative until the 
center of the microburst, i.e., x = 0. Then it is measured positive. To perform 
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Figure 6. Description of coordinate system and approach path orientation 










Plane #1 



Figure 7. Description of coordinate system and data planes 
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a simulation, the runway can be moved back and forth along the path and either 
takeoffs or landings can be simulated at any position relative to the center of 
the microburst 


The important parameters, in terms of assuring that.spatial wind shear is 
incorporated into the aircraft equations of motion, are a and g and angular 
rotation. If the equations are written in the inertial coordinate system, then 
spatial derivatives appear in the equations of motion only in the aerodynamics 
coefficients (Figure 9), in the a and 3 values, and in the angular rotations. 
Consider the a and § terms.as shown in Figure 10. Normally the derivative of 
alpha is taken to provide a as shown in the figure. The values u and w are 
relative velocity values. Beta is handled similarly. To obtain the derivatives 
of the velocity, you write the equations in the^body coordinate system and. 
solve for u and w (see Figure 10). Notice^that u and w contain spatial deriv- 
atives of the wind. Thus when computing a, the spatial derivatives automatically 
appear in your equations of motion. Remember that W^, Wy, and W^ are given in 
earth coordinates and must be transformed to the body coordinates. 


The spatial wind derivatives become even more important in their effect on 
rotations which appear in the aerodynamic coefficients (see Figure 9). For 
example, the lift coefficients, as aerodynamic derivatives with respect to p and 
q, are with respect to relative rotation. What is typically done in most systems 
of equations of motion is to solve for inertial values. Moment equations are 
typically formulated in an inertial coordinate system. What happens, however, is 
that the airplane is rotating relative to the earth (the normally computed rota- 
tion values), but the wind is also rotating (see Figure 11). The JAWS data show 

there is a lot of rotation in the wind (see Figure 12); therefore, the rotation 

of the wind must be subtracted from the inertial rotation of the airplane to get 
the relative rotation. The rotation of the wind comes directly from the nine 
derivatives, because angular rotation of the atmosphere is the anti -symmetric 
part of the tensor gradient (see Figure 13). The equations for the rotation of 
the aircraft relative to the airmass are given in Figure 13. The inertial values 
calculated from the moment equation are reduced by the angular rotation of the 
wind and appropriately transformed into the body coordinate system. These resul- 
ting values should be used to calculate the aerodynamic coefficients. We believe 

that these terms are fairly significant. 

Details of the wind shear models are given in reference 5 and of the aircraft 
equations of motion in reference 6. Standardization of the procedure for imple- 
menting wind shear into simulators is imperative if meaningful training and 
exchange of results are to be achieved. 
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Lift Force: 


L = C L oV 2 A/2 

Ci - Ci + Ci (°0 + Ci (^c) + ( Ci (q) + Ci .(cO) 

L o L a L i E b - v L q L a 

+ ground effect 

Drag Force: 

D = C d pV 2 A/2 

C D = C D 0 + c D a (a * + Cp R ( e R> + c Dg (s) + ground effect 
where 

Cn = drag coefficient due to rudder 
U R 

Cjj = drag coefficient due to side-slip 

B 

Moment: 

M = C m pV 2 C/2 

Cf/j = + Cjv| ( a ) + C M . (5^) + 2 V^C^ (a) + C ri .(«)) 

+ Cp| * e p) + Cf.] (e) + engine thrust effect 
R 8 

+ ground effect 

where 

C Mr = change pitch due to rudder 

C M = change pitch due to side-slip 
"8 


Figure 9. Lift, drag, and moment coefficients 
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Figure 12. Shows rotation of wind field from JAWS data 




H Relative angular rotation vector: 
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Figure 13. Angular rotation relative to atmosphere 
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